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Abstract 


With progressive miniaturization in integration technology in the electronic industry, 
intense efforts are being made in the design and development of printed circuit boards to 
enhance the capability of dissipating the enormous amount of heat flux generated by 
componenis mounted over them. 

The present experimental work has emphatically targeted to simulate printed 
circuit boaro type heated plate configurations. Such configurations are often used in 
numerous eit;ctronic devices like main frame computers, personal computers and more 
sophisticated modern data processing devices. This study has been done with uniform 
surface temperature condition for the heated plate. 

Th ; experimental data which has been acquired in considerable volume helps in 
visualizing the trends in various thermal parameters like surface heat flux, total heat 
dissipated during the cooling, Nusselt number, etc. with relative variation in physical 
parameter.! line plate spacing, channel height etc. These data can be used as a preliminary 
basis for vStimating the cooling performance of a printed circuit board and also for 
designing such boards with given number of heat generating chips and other elements 
over then*,. I’he experimental results have been compared with the theoretical and 
numerical lefjlts for the similar configurations and a few empirical relations for the 
design purposes are suggested. 
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Nomenclature 


At Activation energy 

Bo,,- Boussinesq number measured along y 

Ea Characteristic activation energy 

Gr Grashoff number 

Kb Boltzman constant 

Ni Nusselt number 

Pcj- Peclet number measured along y 

Re Reynolds number 

Ra Rayleigh number 

Q Heat Flux 

T Temperature, K 

X Plate Spacing, mm 

Xop, Optimum plate spacing, mm 

Qop Optimum heat flux, Wcm'^ 

Gravitational acceleration, ms' 


g 



Figure l.Ia 
Figure 1.1 b 
Figure l.I c 
Figure 1.2« 
Figure l.?b 
Figure 1.3 
Figure 2.1 
Figure 2.2 
Figure 2.3 
Figure 2.4 
Figure 2.5 
Figure 3.1 
Figure 3.2 
Figure 3.3 
Figure 3.4 
Figure 3.5 
Figure 3.(i 
Figure 3." 
Figure 3M 
Figure 4.1 
Figure 4..? 



Electronic enclosure 

Chips mounted on the board 

Heat generating chip 

Heat transfer in an electronic component 

Path of Heat transfer 

Flow chart for cooling of a PCB 

Schematic of the test rig used by Tiwari and Jaluria 

Ligranie’s setup for mixed convection in straight and curved channels 

Schematic of the test rig used by Maughan 

Schematic of previous setup used by Maughan 

Schematic of setup used by Ramchandran 

Experimental setup 

Experimental test rig 

Air intake unit 

Test section 

Diffuser 

Suction chamber 
Test section (picture) 

Thermocouple arrangement for heated plate 
Nusselt No. Versus Rayleigh No. 

Heat Flux versus Plate Spacing 



Figure 4.3 
Figure 4.4 
Figure 4.5 
Figure 4.6 
Figure 4.7 
Figure 4.{! 

Figure 4.V 
Figure 4. 10 


Nusselt No. Versus Plate Spacing for different static Pressures (Ts =50°C) 
Nusselt No. Versus Plate Spacing for different static Pressures (Ts =60°C) 
Nusselt No. Versus Plate Spacing for different static Pressures (Ts =70°C) 
Heat Flux Versus Plate Spacing for different static Pressures (Ts =50°C) 

Heat Flux versus Plate Spacing for different static Pressures (Ts =60”C) 

Heat Flux Versus Plate Spacing for different static Pressures (Ts =70°C) 

Optimum Spacing versus Static Pressures for different Plate Surface Temps. 
Optimum Heat Flux versus Static Pressures for different Plate Surface Temps. 



Chapterl 


GENERAL BACKGROUND 


1.1 Introduction 

Electronic cooling problems have always been formidable for thermal engineers with 
respect to their solutions and verifications. This difficulty can be attributed to the 
presence of composite material with different thermal conductivities and multiple 
heat sources with different intensities. 

In electronic cooling the objective is to ensure that component junction 
temperature meets the design specification (80°-100° C) and cooling system is 
effective enough to provide error free operations for worst possible environment. In 
other words, electronic cooling is meant for ensuring reliable operation under 
temperature change. Electronic cooling application areas include PCB’s, IC packages, 
T.V.sets, Digital multimeters etc. 

Miniat'-irization of electronic components has resulted in complex circuits and 
highly dense electronic boards characterized by high rate of heat dissipation per unit 
of comporyent area. Heat flux in modern electronic devices can be of the order of 
about 100 W/cm^. This requires an efficient cooling system with high heat dissipative 
capability. Figures l.la,b&c show atypical electronic enclosure with its components. 





Figurel. Ic; Heat generating chip 


1.2 Effects of Temperature Rise and Failure Modes 

It is important to know why temperature and its reduction are important in electronic 
systems. Temperature induced stresses are of common knowledge. This could be very 
well descri)»ed by fuse burnout case. The fuse bulb can not dissipate the excess heat 
brought about by the power change. Here thermally induced stresses cause the 
mechanical failure, as there is no scope for expansion. Thermally induced stresses 
occur m electronic components when they can not be cooled adequately, resulting in 
eventuai failure. From operational point of view, temperature plays an important role 
(sensitivity of silicon based devices to temperature is a known phenomenon). The 
switching voltage response is highly temperature dependent. Large temperature 
variation at chip level will increase the range over which switching occurs leading to 
switching errors. 

Tonperature induced failures are classified as follows: 

(a) Immediate destruction feilures. 

(b) Increase in the rate of random failures one time(long term 
reliability) 

(c) Operational failure at high temperature. 

(d) Intermittent operational failures at high temperature (parts pf the 
circuit do not work together). 

deduction of the ten^erature at the device level also reduces stress in the 
compo .Menl internal structure. Ten]5)erature also plays a significant role in other feilure 
mechanisms that govern the reliability of the component. This can be very well 



depicted looking at the activation energy and its relation to temperature as 
described by the Arrhenius equation can very well depict this: 

A^=EXP 

Above equation shows an exponential dependence of activation energy 
(signifying failure rate) on temperature. 

1.3 Heat Transfer in Printed Circuit Boards 

The thermal transport is described by coupling between components and system parts. 
To understand electronic cooling, it is important to know about the generation of heat, 
spreading >ind eventual departure at component level. T h is gives us an idea as to how 
the heat gets transferred from the component to the fluid which is an eventual sink for 
the fluid (gas or liquid) cooled system. 

The generated heat at the chip gets transported through any available path 
with th? bulk being transferred through the least resistant one. The sink where heat is 
eventu'illy transferred is cooling fluid. The paths available for heat flow are through 
moldin;;:, rraterial and the leads. Each material regardless of its thickness impedes the 
flow of heat as it travels from sources to the sinks. The wire bond also provides 
another a^'enue for transport of heat. This can become a major path as the cross 
sectional area is proportionally large (typical wire bond diameter is 5|im). 

As soon as heat reaches the leads, part of it is conducted to the board as rest is 
either convected or radiated to the ambient. Combination of multiple heat sources and 
different possible paths for heat flow creates complex and non-uniform temperature 


1 

fi 

n 


1 

U" 

'tJ 




fields. Figures 1.2a and 1.2b show various heat transfers mechanisms and their paths 
followett in a typical PCB. 

1.4 ElficVronic Cooling Methods 

Typical eleclionic cooling methods adopted for heat transfer are as follow: 

(a) Natural Convection Cooling 

Natural or free convection plays a very important role in cooling of electronic 
equipments. In this type of convection cooling, temperature difference in the fluid 
produce density gradients. As a result buoyancy effects drive movement of fluid. 

Natui'a! convection cooling is often used in situations where there is low power 
dissipation or higher temperatures can be allowed. Normally heat transfer coefficient 
achievable in natural convection is of the order of 5 - 25 WW K. This limits the heat 
flux to a. value of 0.1 - 1.0 W/cm^ and allowable temperature to less than 100°C. 
Despite significantly lower surface heat transfer coefficient natural convection is 
preferred tor low end equipments, corq)ared to forced convection or immersion 
boiling, besjause of its inherent reliability and simplicity. 

Blectronic test equipments, consumer electronics, low-end computer packages 
and communication switches are often cooled by natural convection. 

(b) F<krc«d Convection cooling 

In forced convection the fluid movement is due to some external action such as a fan 
or pump. iTiis type of cooling method is used for h^her power dissipation levels or 
in applications such as the space shuttle where near zero gravity eliminates buoyancy 






driven fiow Experiments show that forced convection can enhance heat transfer by 
10 times as compared to natural convection. High heat dissipation leads to increase in 
costs, noise, power consumption and complexity. Forced convection has lead to 
signific;; nt .idvances in the performance of indirect cooling schemes such as micro 
channei coolmg. Typical micro channel with air results in thermal resistance as low as 
OJS^C/W, which can be achieved for airflow velocity of 37m/s. 

forced convection cooling is commonly used for electronic systems with 
high-speed memory chips, multi chip modules etc and other large heat generating 
circuit ebments. 

(c) Mixed Convection Cooling 

Mixed convection is dejSned as a heat transfer situation where both natural convection 
and forced convection heat transfer mechanisms interact. In a vertical passage, the 
internal mifin flow can be characterized in two ways: 

(1) Buoyancy assisted (Upward Flow) 

(2) Buoyancy opposed (Downward Flow) 

In buoyancy assisted flow, natural convection created by buoyancy is m the same 
directicn as the bulk flow. In contrast, the downward flow is based on its direction 
opposho £0 the natural convection. Regardless of the type of convection, the heat 
transfer is always between the vertical plate and the fluid reservoir between them. 
Non dimensional parameter that determmes the type of convection dominating is 
given 03 fcilows: 
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In past few years mixed convection in a heated channel has received 
considerable attention due to its extensive practical applications, such applications 
include house hold electronic devices and equipments, industrial electronic packages, 
large mainframe systems, avionics and space electronics, nuclear reactor electronic 
components and so on. Here in all cases, the target is to optimize the power 
consumption and heat flux dissipation as well. It is found that enhancement of heat 
transfer is most pronounced when the channel is horizontal and heat transfer is taking 
place from the bottom surface, yet in electronic devices vertical plate stacking is 
preferred because of space limitations to allow heat to be conducted from both sides. 
1.5 Methodology for designing a cooling system 

From simplicity point of view, we consider a single chip module residing on a board 
and air is flowing over the component for the pxirpose of maintaining its junction 
temperature below the design specification. The objective here is to obtain the 
junction temperature of the component, which is a function of the board thermal 
conductivity, component surfece and fluid temperatures. The methodology adopted 
here is shown in the following figurel.3. 
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Chapter 2 


LITERATURE SURVEY 


2.1 Introduction 

This chapter is a brief account of works done in the area of electronic cooling. This 
includes the summary of analytical, numerical and experimental works reported in 
literature. This literature survey is not only to fecilitate the improvement of existing 
scheme:^ but also to provide a basis for development of new ones. The survey is not 
only linuted to the investigation performed over the last few years, but it also covers 
some cf the landmark works done over the few decades as well. It includes 
informa tior. on possible simulation of fluid flow and heat transfer processes 
specifically related to the thermal control of the semiconductor devices, modules and 
the cor.'*plele systems. 

Various research efforts done in this field comprise of the following heat transfer 




mechanisms. 


1. Natural Convection. 

2. Forced Convection. 

3. Mixed Convection. 

4. Phase change. 

5. Jet Impingement. 

• Each of the mechanisms have been explored experimentally, analytically and 
also numerically. 

This chapter deals with mixed convection only and therefore works done in this area 
only are summarized. 

2.2 ANALYTICAL STUDIES 

Min-Joon et al. [1] investigated one dimensional mixed convection of liquid metal with 
forced circulation and unequal heating. The steady solutions of the corresponding heat 
equations with proper boundry conditions were obtained analytically. The solutions 
showed formation of local circulation loop by unequal heating and deformation of the 
circulation loop by the forced convection. It was also found that the heat transfer 
characteristics were almost independent of the unequal heating conditions although 
the deformation occured more easily for a more symmetric case. 

2.3 EXPERIMENTAL STUDIES 

Tewari and Jaluria [2] carried out an experimental study on the fundamental aspects of 



the conjugate, mixed convective heat transfer from the finite width heat sources which 
were of negligible thickness and had uniform heat flux input at the surface. These were 
located on a horizontally or vertically oriented fiat plate. Two dimensional flow cir- 
cumstances were simulated and the convective heat transfer mechanism was studied 
in detail by measuring surface temperature and determining the heat transfer coef- 
ficient for two heated strips representing isothermal sources. Results obtained have 
great relevance in cooling of electronic systems. Results showed sharp variations in the 
surface temperatures near the leading and trailing edges of the heated stripes for both 
horizontal and vertical surface orientations. It was also found that in the horizontal 
orientation of the surface even a strong heat source has a negligible effect on the sur- 
face temperature of relatively weak neighbouring heat source. This was attributed to 
the insignificant interaction between vertical rising wakes- from the two heat sources. 
Schematic of the setup used is shown in figure 2.1. 

Chen et al. [3] investigated the linear stability analysis of mixed convection in a differ- 
entially heated vertical channel for Prandtl number values of 0.7, 7.0, 100, and 1000. 
The viscous flow investigated in this paper was the mixed convection, which was driven 
by an external pressure gradient and also by a buoyancy force, between the two par- 
allel long vertical plates separated by a distance. The gravitational force was aligned 
in the negative x-direction. There was a fixed temperature difference, produced by 
maintaining the two vertical walls at different temperatures T1 and T2 respectively. 
It was found that both the Prandtl number and the Reynolds number were important 
parameters in determining the critical Grashoff number Grc, critical wave number ckc, 
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Figure 2.1; Schematic of the test rig used by Tiwari and Jaluria 





wave speed Cr and instability mechanism for higher Prandtl numbers. For the lower 
Prandtl numbers, the effects of the Prandtl number and the Reynolds number were 
relatively small. The most significant finding was that the local minimum wave num- 
ber was 8.0 for Pr=1000.0, which was substantially higher than those found before 
for earlier heated flows. In general, for mixed convection channel flows, the instability 
characteristics of the differentially heated flows were found to be substantially differ- 
ent from those of uniformly heated flows. For constant heat flux case, instability type 
was in general Thermal buoyant instability for the buoyancy assisted flow. In case of 
constant wall temperature, thermal shear instability was generally found for low to 
moderate Reynolds numbers. 

Tso et al. [4] measured convective heat transfer from a linear array of flush mounted 
heat source in a vertical flow up channel. Water was used to simulate the dielectric 
fluid. The dimensions of the simulated chips were held constant while height of the 
channel was varied. The experiments were conducted in a closed-loop liquid cooling 
flow facility, with a vertical up-flow, in a plexiglas test section. The first chip was 
located 680.0 mm downstream of the channel inlet, providing a minimum hydrody- 
namic entry length of 50 hydraulic diameters. Channel height was varied over values 
of 0.5, 0.7, 1.0 times the heated source length. The heat flux was set at three values 
of 0.5, 10.0 and 20.0 W/crn^. Reynolds number based on heat source length ranged 
from 6 X 10^ to 8 x 10^. The results obtained agree well with the results for air-cooling 
of heated sources of similar geometry when Peclet numbers were used to correlate 
the data in the laminar flow and the Nusselt numbers were normalized against the 



Prandtl numbers in the turbulent flow. They suggested that data from air cooling 
might be used to predict the heat transfer characteristics of the liquid cooling of a 
similar geometry if the Prandtl number scaling was considered. 

Tseng et al. [5] carried out an experimental investigation combining flow visualisation 
and temperature requirements to unravel the effect of tapering the top plate of the 
rectangular duct about the detailed changes in the vortex flow structure in a mixed 
convection air flow through a bottom heated trapezoidal duct. Special care was taken 
on the flow stabilization and the weakening of the vortex flow by the acceleration of the 
primary flow for the various Reynolds and Grashof numbers. The spanwise-averaged 
Nusselt numbers for the horizontal rectangular and tapering ducts were obtained to 
further detect the temporal stabiltiy of the flow. Over the ranges of the Re and Gr 
investigated for 5 < i?e < 102 and 1.0 x 10'^ <Gr< 1.7 x 10®, the vortex flow induced 
in the rectangular duct exhibited temporal transition from a steady laminar to the 
time periodic and then to chaotic state at incerasing buoyancy to inertia ratio. During 
the course of the study it was realized that although the top plate tapering resulted in 
the complete stabilization of the vortex flow, the weaker vortices still dominated the 
secondary flow in the tapering duct. 

Zhang et al. [6] conducted experiments and analysis of buoyancy-assisted mixed con- 
vection in a vertical square channel with asymmetric heating conditions. Most pub- 
lished literatures in mixed convection are based on test geometry with all sides heated 
. This work was an assisted mixed convection in a vertical square channel under 
four symmetric heating conditions i.e. with one or two sides heated. This created a 
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stronger three dimensional flow disturbance that could change heat flow characteris- 
tics both upstream and downstream. The paper covers the laminar, trasition, fully 
turblent flow regimes. The Reynolds number was varied from 200 to 11200 and the 
buoyancy parameter, was varied from 0.02 to 200. It was found that local heat 
transfer coefficient increased with the increasing buoyancy parameter. The Nusselt 
number ratio, was observed to increase significantly with the increased turbulent 
mixing caused by more irregular heating conditions. The increase in the coefficient 
could be attributed to three-dimensional large scale buoyancy- induced motion in the 
flow field. 

Ligrani et al. [7] investigated mixed convection in straight and curved channels when 
buoyant forces were perpendicular to the bulk flow direction. Grashof number based 
on the temperature difference ranged from 0 to 50 x 10®, and the Reynolds number 
based on channel width ranged from 343 to 2057 which gave laminar and transition 
forced flow. The flow behaviour was classified into four different regimes viz, 

1. pure natural convection, 

2. mixed convection with strong buoyancy, 

3. mixed convection with weak buoyancy, and 

4. pure forced convection. 

The authors discussed the difference in the behaviour between these different regimes 
along with the curvature influences and suggested Nusselt number correlations for the 
two mixed convection regimes. Test rig used for the purpose is shown in figure 2.2. 
Gau et al. [8] conducted an experimental study of the mixed covection in a heated 
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divergent channel with one wall held in the vertical postion. The vertical wall was 
heated uniformly while, the opposite wall, which had a divergent angle of 3.0 degree 
was tilted away from heated wall and was well insulated. They studied the occurence 
and the structure of the reversed flow and the heat transfer in a heated vertical di- 
vergent channel. For the purpose of comparison, results for the parallel plate channel 
were also presented and discussed. The Reynolds number for the main forced flow 
ranged from 100 to 4000 and the buoyancy parameter, varied from 0.3 to 907. 
Flow reversal was found for assisted and opposed convection. It was found that the 
divergence of the channel decelerated the main stream such that flow reversed at a 
much lower buoyancy parameter. The adverse pressure gradient tended to push the 
reversed flow upstream and led to a deeper penetration of the reversed flow into the 
channel. The destablization effect of the divergent channel could lead to breakdown 
of vortices and transition to turbulent flow. This could significantly enhance the heat 
transfer. Temperature fluctuation measurements at different locations were used to 
indicate oscillations and fluctuations of the reversed flow. The effect of the buoyancy 
parameter on the Nusselt number and the reversed flow structure was discussed. It was 
found that for the buoyancy assisted convection Nusselt number increased with the 
buoyancy parameter. For opposed convection the increas in Nu with the buoyancy 
parameter was related to the occurrence of reversed flow which interacted strongly 
with the mainstream. The average Nusselt number was determined and correlated in 
terms of relevant nondimensional parameters for pure, forced and mixed convection 
respectively. The fluid velocity in the channel was measured with a TSI hot-wire 
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anemometer using an improved calibration system. 

Maughan et al. [9] investigated heat transfer enhancement in a mixed convection flow 
between parallel plates heated uniformly from below. They performed both the ex- 
perimental and numerical study. Flow visualization and heat transfer measurements 
for air flow in a horizontal channel were performed over a wide range of conditions to 
delineate the regions of forced and mixed convection. The onset of secondary flow was 
found to precede appreciable heat transfer enhancement. Except for large values of 
Nusselt number prior to the onset of enhancement could be predicted with forced 
convection correlations. Their numerical computations showed heat transfer enhance- 
ment prior to the onset of thermal instability in a horizontal channel due to fluid 
expansion and an induced pressure gradient, and in an inclined channel due to com- 
ponent of gravity in the flow direction. In the numerical work, the plate conduction 
eflfects were neglected, and the thermal entry region was investigated by prescribing 
laminar fully developed flow at the start of the heated section. Assuming steady two 
dimensional laminar flow and constant properties and the Bousinesq approximation, 
the appropriate governing equations were non-dimensionalized. The experiments were 
performed using the horizontal air channel which along with other accessories is shown 
in figure 2.3. 

Ambient air was drawn into blower, metered and conditioned before passing through 
the channel. Warm air leaving the blower was cooled to ambient temperature with a 
compact air to water heat exchanger. The inlet temperature of water was controlled 
by a data acquisition computer which based upon the inputs provided by the ther- 
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Figure 2.3: Schematic of test rig used by Maughan 
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mocouples in the air stream leaving the exchanger, opened or closed an electronically 
controlled valve mixing hot and cold water. A uniform and constant inlet temperature 
was assured with additional heat exchanger located just prior to the channel entrance. 
This exchanger was driven by a constant temperature bath held at the ambient tem- 
perature. 

Maughan et al. [10] performed few experiments to investigate mixed convection heat 
transfer in the thermal entry region of a parallel plate channel heated uniformly from 
below. The effect of surface heat flux and the channel orientation on the local Nus- 
selt number was studied for Pr=0.7, 125 < Re < 500, 7 x 10^ < Gr* < 10®, and 
0 < 9 < 30°. Heat transfer initially was dominated by forced convection and showed 
a rapid decline in the Nusselt number. Following the onset of thermal instability, 
secondary flow development caused a sharp increase in the Nusselt number, which 
was followed by maximum subsequent oscillations. The Nusselt number oscillations 
eventually decayed, yielding a fullly developed value which depended on the Grashof 
number. The onset of instability was delayed by decreasing the Grashof number and/or 
by increasing the Reynolds number and the inclination angle. For the inclined channel, 
significant heat transfer enhancement occurred prior to onset of the secondary flow. 
The experimental test rig used is shown in the figure 2.4. 

Ramchandran et al. [11] conducted measurements and prediction of laminar mixed, 
forced and free convection air flow adjacent to an isothermally heated vertical flat sur- 
face. Their experimental test rig is shown in figure 2.5. It basically consisted of low 
turbulence open circuit wind tunnel which could be rotated and fixed at any desired 
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Figure 2.4: Schematic of previous setup used by Maughan 
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Figure 2.5: Schematic of setup used by Ramchandran to study laminar mixed, forced 
and free convection 
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inclination angle. The tunnel had a smooth converging nozzle with a contraction angle 
ratio of 9:1, a straight test section and a smooth diverging diffuser. Plastic honeycomb 
material ( 10.00 cm thick ) and several wire screens were used in the front section of 
the tunnel and along the nozzle length and diffuser to straighten the flow and reduce 
the turbulence level in the test section. The heated flat plate consisted four layers 
which were held together by screws and instruments to provide an isothermal heated 
test surface. In this study local Nusselt numbers and the velocity and temperature 
distributions were presented for both the buoyancy assisting and opposing flow cases 
over the entire mixed convection regime, from pure forced convection limit ( buoyancy 
parameter ^ ^ = 0 ) to the pure free convection limit = inf) . It was found that 

the mixed convection Nusselt numbers were larger than the corresponding pure forced 
and pure free convection limits for buoyancy assisted flow and were smaller than these 
limits for the opposing flow. Flow reversal near the wall was detected for buoyancy 
opposing flow case at a buoyancy parameter of about ^ = 0.20. 

Yadav et al.[12] conducted an experiment to determine heat transfer characteristics 
of a stack of vertical heated plates. The geometrical configuration designed and fab- 
ricated for this work simulated the printed circuit boards(PCB’s) used in various 
electronic devices. The uniform heat flux model was adopted. Heat flux from the 
plates was kept every time at a value such that temperature at any point shall not 
exceed the preferred temperature limit. It was shown that heat transfer rate could 
not be enhanced indefinitely by increasing the coolant driving pressure Certain points 
were identified on the plate surfaces where high heat flux could be safely reached. 
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Several other noticeable trends were observed regarding power consumption etc. 

2.4 NUMERICAL STUDIES 

Rao ot al.[13] numerically studied a two dimensional steady state conjugate problem 
with laminar compressible flow in mixed convection with surface radiation from a ver- 
tical plate with a flush mounted discrete heat source. The effects of the magnitude 
and location of the heat source, the material and surface properties of the plate and 
the free stream velocity on both heat transfer and fluid flow were studied. Vorticity 
stream function was used to write the governing equations and were finally solved 
using finite volume based difference method. It was found out that the best position 
for the heat source was the leading edge of the plate and the least preferable position 
was the trailing edge, as the latter position increased 0(non dimensional temperature) 
by 35 - 50 percent for the entire range of Rn. For a given Rn =1, convection was 
dominant taking away as much as 90 - 95 percent of heat generated for good reflecting 
surfaces. 

Papanicolaou and Jaluria [14] conducted a numerical study on the combined forced 
and natural convection air cooling of heat dissipating electronic components located 
in a rectangular enclosure. Flow and thermal fields in flowing air, at the walls of 
the enclosure and mounted electronic components were obtained. Results for laminar 
and oscillatory domains were presented which have vital importance in the search for 
a suitable placement of electronic components in a given enclosure for effective heat 
removal. 

Davalath and Bayazetogue [15] investigated conjugate heat transfer for two dimen- 
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sional developing flow over an array of rectangular blocks which represented finite 
heat sources on parallel plates. Incompressible flow over multiple blocks was modeled 
using fully elliptic forms of Navier Stokes equations. A control volume based finite 
difference procedure with appropriate averaging for diffusion coefficient was used to 
handle coupling between solid and fluid regions. They have extended the analysis to 
study the optimum spacing for a fixed heat input to achieve desired maximum tem- 
perature at the heat source. 

Barletta et al. [16] investigated mixed convection of the power-law fluid in a verti- 
cal annular duct in a regime of laminar fully developed flow. Uniform and unequal 
temperature were prescribed on the inner and outer boundary walls. The momentum 
balance and the energy balance equations as well as the viscous stress constitutive 
equations were solved analytically in order to obtain the velocity fields, the viscous 
stress field and the temperature field. The analysis of the cases was done for : 


1 . 

2 . 


Mixed convection of a Newtonian fluid, and 


Forced convection of the power law fluid given by 

dvx |n-l rfv® 


Tyx — ^1 cly I dy 


for n = 1, this equation stands for Newtonian Fluids. 


In the first case it was shown that the condition for the onset of the reversed flow 
implied the existence of a threshold value for A =| | which depended on the 7 = || 

(ratio of inner to outer radii). When These Threshold values were exceeded, flow 
reversal occurred next to the inner boundary wall. In the case of mixed convection 
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for the power law fluids, the flow reversal conditions were investigated. As for the 
Newtonian fluids, flow reversal occui'red when the threshold values of modA were ex- 
ceeded. These threshold values depended on both 7 and inverse of power law index 
(m). Further for the fixed values of 7 , Threshold values were increasing function of m. 
Moukallad et al. [17] numerically studied the mixed convection heat transfer in a 
channel with heated curved surface bounded by a vertical adiabatic wall. Two cases 
were considered, in the first case the flow experienced a convex curvature on an in- 
creasing cross-sectional flow area ( adverse pressure gradient ), while in the second 
case, the flow experienced a concave curvature with a decreasing flow cross-section ( 
favorable pressure gradient ). Governing equations were conservation of mass, mo- 
mentum and energy. Flow was considered to be steady, laminar and two dimensional. 
The volume approach was adopted to solve numerically the coupled system of equa- 
tions governing the flow and temperature fields. For the convex entry channel, it was 
observed that for low [Gr/Re^), the Nusselt number reached a minimum near the 
onset of flow seperation along the heated wall of the convex-entry channel and then 
decreased towards the channel exit. The effect of increasing Prandtl number on the 
total heat transfer was considerably greater than the effect of increasing (Gr/Re^). 
For concave-entry channels, the flow seperation was not encountered. The decreasing 
cross-sectional area combined with buoyancy increased the channel velocities near the 
wall gradients. Therefore, the surface heat transfer rates were considerably greater 
than in the convex-entry channel. The overall heat transfer in a concave-entry chan- 
nel was always greater than a straight channel of equal heated surface area, a critical 
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vah'.e :or the ratio (Gr/Re^) existed, below which heat transfer enhancement was ob- 
tain ed with concave-entry channels. 

Yadav and Kant [18] developed a theoretical model to numerically simulate mixed 
convection cooling in a stack of vertical heated plates representing the PCBs. For 
the gi' an spacing between the plates, they estimated the maximum amount of heat 
that covld possibly be transfered to the coolant at a certain static pressure difference 
between the top and the bottom of the stack and also the maximum possible heat flux 
from the surface being cooled. It was shown that the need for the optimum spacing 
betwe<}n the plates became more important when the driving pressure reached a rel- 
atively high level. The following empirical relation was suggested for predicting the 
heat iiux corresponding to the optimum plate spacing for a given static pressure : 

Qi = ‘oa^ + ^ - 0-205 

Chong et al. [19] looked upon the mixed convection in a differentially heated vertical 
channel for various Prandtl numbers. Galerkin method was used to solve coupled mo- 
mej'.tum and energy, equations along with continuity equation. They predicted that 
Prandil and Reynolds numbers had very important effect on the Grashof number, 
wa'\ 3 number, wave speed, instability mechanism for higher Prandtl numbers. It was 
found 'hat the local minimum wave numbers were as high as 8 for Pr=100, which was 
substantially higher than those found earlier for other heated flows. 

2.5 CLOSURE 

Research efforts in the area of cooling of electronic components have used several heat 
transfer mechanisms as mentioned in the text. Mixed convection happens to be most 
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common and convinient due to certain merits. Due to the miniaturizatuion of elec- 
tronic components, heat dissipation rates per unit area are quite large and the efforts 
are being mad(^ to increase it even more by using one or a combination of mechanisms. 
Lot of work has been done in the field of mixed convection as such. Experimental 
studies have been made with staright, curved and divergent channels using constant 
heat fills and constant surface temperature conditions. Buoyancy and pressure driven 
flows have been studied. But the subdomain of these works relating to the mixed 
convection cooling in electronic equipments is not very large. The models adopted 
for these specific studies have been relooked and modified according the recent needs. 
Selected works pertaining to the electronic cooling were presented as briefly as possi- 
ble without distorting the essence of the original work. Going through the presented 
material one can realize the need of target oriented reaserch work ie., the modelling 
of printed circuit boards and their possible configurations inside the cabinet etc. 
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Chapter 3 


PROBLEM DESCRIPTION AND 
EXPERIMENTAL SET-UP 


3.1 Introduction 

As already discussed, mixed convection cooling has a clear edge over natural 
convection and forced convection cooling methods. It possesses the simplicity and 
less noisy features of natural convection simultaneously having high heat removal 
capability of forced convection. In order to simulate a mixed convection cooling 
problem one has to take into consideration various parameters that go into making of 
an actual heat generating source in today’s PCB’s. Some of them are like the 
alignment of heating element, its location with respect to the cooling system and 
cooling fluid, type of cooling fluid used and so on. For such a complex problem like 
the present one, it becomes necessary to follow a systematic approach and make 
assumptions wherever necessary. 




3.2 Experimental Requirements 

In order to conduct experimental study on the cooling mechanisms inside various 
electronic equipments and devices, one has to model a PCB along the following lines. 
As already discussed, there are various parameters which should be taken into 
account while simulating a true printed circuit board. 

(1) It’s possible arrangement. 

(2) It’s possible variety of dimensions. 

(3) Possible orientation of a single PCB with respect to other PCBs. 

(4) Variety of possible mounting or arrangements for electronic chips having 
different temperatures. 

(5) Type of cooling fluid used (air or water). 

(6) Modes of cooling, whether free or forced or combination of both. 

(7) The unit driving the working fluid or static pressure creating device and its 
location. 

But in the current experimental study it is not possible to have aU the 
flexibility of parameters. We shall stick to certain parameters, which satisfy our 
requirements for obtaining the maximum cooling effects and least power 
consumption. Looking at the current requirements of electronic cooling following 
parameters are taken care of: 

(1) Material used must possess the specified thermodynamic properties within 


the specified temperature range. 



(2) Driving fan must be located at a place from where the acoustic 
electromagnetic disturbances shall have minimum effect, as it can affect 
badly the sensitive components mounted inside. 

(3) Working fluid must move within the specified limits as turbulence may 
lead to damage. 

(4) There should be provision for air filters to remove unwanted dust particles, 
which could prove to be damaging agents for delicate components and 
may choke the micro passages for coolant entry. 

3.3 Objectives 

The objectives set for the present study were as follows: 

(1 ) Design and development of an experimental apparatus for Simulation of Mixed 
Convection Cooling of PCB type configurations. 

(2) Make a parametric study of the variables involved e.g. inter plate spacing and 
static pressure to determine Convective heat transfer coefficient, average and 
maximum heat flux. 

(3) Suggest suitable correlation among the variables. 

3.4 Experimental Set-up 

In the present study mixed convection cooling was simulated eiqperimentally for PCB 
type configurations. In this experiment one vertical channel formed by two simulated 
PCBs was placed in the test section where a lead screw mechanism vras used for 
varying the spacing. A suitable static pressure creating device (a fen unit) was 
provided to suck the air passing through the vertical channel- One of the two plates 
was acting as an insulating wall while the other consisted of heating elements. There 



were five heating elements in the heated plate. Vertical orientation of the channel 
helped in achieving the mixed convection. Since the suction device was placed at the 
top of the test section, the buoyancy assisted mixed convection flow was achieved. 
Completely assembled experimental test setup is shown in figure 3.1. The objectives 
of the experiment dictated the test rig configuration. It was essential to take into 
account the accessibility and the flexibility of its components so that there was not 
much difficulty in taking care of maintenance and repair problems. Provisions for 
rapid configurational changes were very important and these were incorporated in the 
design. Figure 3.2 shows the schematic diagram of the test rig designed for the 
purpose. 

Test rig consisted of the following major parts: 

(a) Air intake unit, 

(b) Test section, 

(c) Diffuser unit, and 

(d) Suction chamber (fen unit). 

The detailed description as follows: 

(a) Air intake unit 

This unit helped in achieving smooth entry of air at the test rig inlet. This unit 
was designed such that the flow at the inlet of test section was fully developed 
as fer as possible. It was basically a converging duct with rectangular cross- 
section. The dimensions are shown in figure 3.3. The material used for intake 
unit was galvanized iron sheet. Since few thermocouples were placed at the 




Figure 3.1: Experimental Setup 






Figure 3.2: Experimental test-rigC 




Figure 3.3: Air intake unit 



exit of this unit, the metallic sheet helped in eliminating the electrical 
disturbances due to other electrical appliances used in the lab. 

(b) Test section 

The cross-sectional area of test section primarily determines overall size of the 
test .section. Te.st section being the most important part was designed taking 
into account structural cost, power cost etc. Many shapes for the cross section 
of test section are used like round, elliptical, square, hexagonal and octagonal 
with minimal difference in losses due to its shape. But the current study used a 
square cross section. Length of the test section should be optimized as it 
affected the power consumption. It should not exceed the required length. 
Figure 3.4 shows the test section with its actual dimensions. Fifteen- 
millimeter thick glass plate was used for making the test section. Glass being a 
good insulator was a right choice for this purpose, it also provided a window 
for viewing the whole plate assembly. 

(c) Diffuser 

Material chosen for the diffuser was perspex. Figure 3.5 shows detailed 
dimensional diffuser. Diffuser was designed in order to avoid separation of 
flow, which would have otherwise led to vibrations, oscillating fan-loading 
variation in tests section velocities. 

(d) Suction Chamber 

In order to create a pressure differential across the channel, suction chamber 
was provided with four small fans with the following rating. 

Discharge Rate; 0.2 m^/s 



Teflon Sheet 
Bakelite Sheet 






Rated Speed; 1440 r.p.m. 

The idea behind providing four small fans was to avoid inherent vibration 
probknts in a single fan of large discharge rate. Figure 3.6 shows the 
dimensional wooden box that was used to house the exhaust fens. 

3.4.1 Heated Plates 

While designing a heated plate for the simulation of a PCB, following considerations 
were made: 

{ 1 ) The material of the plate was electrically insulating. 

(2) It had the ability to withstand high temperatures (200 -300°C). 

(3) It should be easily machinable. 

Looking at the above requirements we have selected Teflon as material 
for the plate. Its some useful properties are shown in the table 3.1. Complete 
picture of the plate assembly inside the test section is shown in the Figure3.7. 
Heating elements with complete details are shown in Figure 3.8. Five such 
heating elements were placed between the plates. 

Teflon Properties 
Table 3.1 


Melting Point 

267 ‘'C 

Upper service Temp 

150 ‘^C 

Tensile Strength 

1 

44.5 MPa 

Thermal Conductivity 

0.2 W/m.K 





Figure 3.7: Test Section 
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All dimensions in mm. 


Figure 3.8: Thermocouple arrangement for heated plate 



3.5 Thermocouple Arrangement 

3.5.1 On the plate 

In'jr nicaMiring temperature on the plate surface, a total of ten theiirocouples 
%scre used, two ibr each section of the plate corresponding to heating 
elements. 

3.5.2 Inside the plate 

Additional ten thermocouples were used in a similar fashion to measure the 
temperature inside the plate assembly. Above temperature measurements were 
u,sed for finding out heat fluxes from the plates. 

In the present work Chromel-alumel thermocouples (K-Type) 
were used as temperature sensors. These thermocouples were used because of 
its availability and large temperature range (-150‘’C - 1350°C). , 

3.6 Experimental Procedure 

The experiments were conducted to obtain heat transfer data for a vea'cal channel 
formed by a printed circuit board (PCB) like vertical configuration and a thermally 
insulated wall. For the modeling of a PCB as a plate with uniform surface temperature, 
different approaches are possible. The approach adopted here can be described as 
follows. 

The entire plate surj^:^ was virtually divided into five zones placed vertically 
one over another. Since the plate height and width chosen were 50 cm and 20 cm 
respectively, each zone ha d dimensions of 10 cm x 20 cm. As the objective was to 
achieve uniform te mp erature over the entire plate Le. same temperature for all the five 
zones, each zone was heated hy a separate stnp-heatmg element. By regulating the 



electrical fxiwer fed to each element, it was possible to maintain a fixed temperature 
over each /tmal surtace and hence ^in a sense) over the entire plate, while it was 
subjected to various cooling rates. As mentioned earlier, four thermocouples were 
nwunted over centerline of each zonal surface, two at the top surface and two beneath it. 
The difference in the temperature given by each pair was used to calculate the heat flux 
value for that zone. Imither, the pair of thermocouples mounted on the strface of the 
zone was used to monitor the temperature of the surface and thus regulating the power 
input to the associated beating element. The uniform surface temperature condition 
clearly meant different values of heat flux for various zones. The experiment was run for 
different value.s of airflow rates (i.e. different static pressures), plate spacinji^s and plate 
surf^e temperatures. A typical experimental run consisted of the following steps: 

1 . The Power was switched on. 

2. The Plate spacing was fixed at a required value. 

3. Mow velocity was adjusted by controlling the power input to the four fljif mounted 
at the top of the apparatus to obtain the required flow rate and the static pressure. 

4. Temperatures on the plate surlEace were monitored using temperature controllers 
and power to the heaters was adjusted until same steady plate surface Lenjperatures 
were obtained in each zone. These steady state temperatures were recorded using 
thermocouples. 


5. Relevant parameters were coirputed. 



Chapter 4 


RESULTS AND DISCUSSION 


4.1 Introduction 

Previous chapter described the test rig and its instrumentation to study heat transfer for 
PCB type configuration. A large volume of data was collecte;! f\»r the uniform 
temperature of the plate simulating the PCB. This chapter is intended to present the 
analysis of the same data. Some remarkable trends were observed :n the heat transfer 
parameters like heat flux, Nusselt number profile etc.; these are beinji highlighted at the 
appropriate places. Empirical relations were also obtained and described here for the 
purpose. 

4.2 Data Range 

Parameter Values for which data were obtained 

1 .Surface Temperature 50, 60 and 70°C 

2.Plate Spacings 15, 20, 30, 40, 50, 60 and '70mm 

1.6, 2.4, 5.4, 9.6 and 15 Pa 


S.Static Pressure Difference 



4.3 Heat Dissipation by Natural Convection 
43.1 Nusselt Number Versus Rayleigh Number 

Figure 4.1 shows variation in Nusseit number with Rayleigh number for pure natural 
convection situation for different uniform surface temperatures. The plot clearly shows a 
substantial increase in Nusseit Number as uniform surface temperature was increased. 
This increas ,* in Nusseit Number was attributed primarily to high heat input required to 
maintain high surface temperatures. 

On ’.he other hand for a given uniform surface temperature there was no 
substantial v ariation in the value of Nusseit Number with respect to Rayleigh Number 
(signifying he e.ffect of buoyancy) except at the lower range of the data. These relatively 
high values of Nusseit Numbers may be due to high-induced flow at lower plate spacings. 
Increase in plate spacing slightly decreased the induced flow velocity, which resulted in a 
decreasing trend of the curve and finally approaching a constant value of the Nusseit 
Number. 

4.3.2 Effect of Plate Spacing 

Figure 4.2 shows the variation of plate heat flux with respect to plate spacing at different 
uniform sui face temperatures. The plot clearly shows the increase in heat flux value as 
surface temperature was increased. This increase in heat flux value was attributed to high 
heat input required in order to maintain high surface temperature. 

As can be seen from the figure, for a given surface temperature, heat flux value 
does not vary much .It shows that after certain spacing heat flux can not be enhanced just 
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Fig 4.1: Nusselt No. versus Rayleigh No. 
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Fig 4.2: Heat Flux versus Plate Spacing 





by increasing the spacing, {'he figure shows that there is a limitation up to which heat can 
be dissipatc<i oy natural coinection. In order to enhance it fiirther one has to resort to 
some other cooling methods. One of them which can help increase the heat dissipation is 
Mixed convection i.e. heat transfer by using both natural as well as forced convection by 
providing suflicient static pressure across the channel to be cooled. 

4.4 Heat Dissipation by Mixed Convection 
4.4.1 E0ect of Plate Spacing on Nusselt Number 

Figure 4.3 ih.iw.s the variation of the Nusselt number with plate spacing for different 
values of tiie :an static pressures. These data were taken for plate surface temperature 
Ts=50°C. l! cf..i be observed jfrom figure 4.3 that for the static pressure values at the 
lower side. Nuiselt number almost remained around 500. As the Static pressure value 
was increa :ed above 5.0 Pa, Nusselt number no longer remained constant for small plate 
spacings ;say up to, 30 mm) beyond which it increased and then decreased and 
approached a constant value as the spacing was increased. It may be noted that for higher 
static pressures, the initial slopes of curves tended to increase. For UWT of 50°C, 
maximum ^^chievable value of Nusselt number was around 5400, which corresponded to 
a static pressure of 1 5 Pa and plate spacing of about 15 mm. 

Figure 4.4 shows the variation of Nusselt number with plate spacing at wall teir^jerature 
of 60°C. The Nusselt number had a slightly higher value at lower spacing and lower static 
pressure atv can be seen by comparing the curve for 2.4 Pa in this plot with that shown in 
figure 4.3. The Nusselt number invariably tended to achieve constant values for various 
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Fin 4.3; Nusselt No. versus Plate Spacing for different Static Pressures (Ts=50^C) 






static pressures. The maximum value of Nusselt number v/as about 6000 for static 
pressure ot 1 .^.0 Pa and plate spacing of 20 mm, however maximum limiting value of 
Nusselt nun'.bor was around 3800. 


figure 4.5 £'.!$(< shows variation of Nusselt number with plate spacing for different static 
pressures at tlie UWT of 70^C. It can be clearly seen that Nusselt number invariably 
increased w .tb the increase in static pressure. As compared to the characteristics shovra in 
figures 4.4 and 4.3, some increase in Nusselt number could be seen even at lower 
pressure valui.s like 1.6 Pa. At smaller spacing, the Nusselt number value was 
substantial! higher, but it approached a lower limiting value as the plate spacing was 
increased. The Maximum values of Nusselt number were noted at plate spacing around 
15 mm icr d’lferent static pressures. For UWT of 70°C, overall maximum value of 
Nusselt nu nbe*’ was around 4400 and the steady state maximum value was about 4000. 

4.4.2 Elfect of Plate Spacing on Heat Flux 

Figure 4.6 she -vs variation in heat flux (W/ern^) with plate spacing and varying static 
pressures for a plate temperature of 50°C. The heat flux undoubtedly increased with an 
increase m static pressure, but the effect of plate spacing variation over heat flux was 
difficult to analyze in a straightforward manner. One thing that was clearly observed was 
that for some initial values of spacing (between 15mm and 40mm), the flux approached 
some peak value before showing oscillations. In the parametric range under consideration 
the maximum value of heat flux occurred at static pressure of 15.0 Pa and spacing of 
around 30 mm. 
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Fig 4.5: Nusselt No. versus Plate Spacing for different Static Pressures (Ts=70°C') 






Figure 4.7 shows the heat flux variation with the plate spacing for different values of 
static pressure I'he UWT was maintained at 60°C. Heat flux invariably increased with 
increase in static pressure. Heat flux values were highest at the minimum value of plate 
spacing (near 15 mm). The heat flux first decreased sharply with increase in plate spacing 
and then attr.ined a constant value. The apparent increase in heat flux value for the larger 
plate spacing ^vas mainly due to some instrumentation and measmement inaccuracies 
rather than due to .some physical reason. The maximum value of heat flux was about 0.35 
W/cm^, which could be achieved for static pressures between 10 and 15 Pa depending on 
the plate sj^acing. For the larger plate spacing the steady state value was near about 0.3 
W/cml 

Figure 4.8 -hews the variation of heat flux with plate spacing for different values of static 
pressure. 1 he (JWT was maintained at 70°C. The pattern followed by the heat flux was 
more or le ss similar to figure 4.7. The value of heat flux invariably increased with the 
increase i.v static pressure. At the low static pressures, the variation in heat flux for 
different plate ..pacing was relatively low and gradually increased with increase in it. The 
maximum value of heat flux achieved at UWT of 70®C was about 0.45 WW at a plate 
spacing of about 15 mm and at a static pressure of 15.0 Pa. The maximum steady state 
value for the heat flux was around 0.43 W/cm^ for the same static pressure. 
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Fig 4.7: Heat Flux versus Plate Spacing for different Static Pressures (Ts=60°C) 






4.5 Effect of Static Pressure on Optimum Spacing 

By optimum spacing we precisely mean that the value of plate spacing for which either 
the value of h- tat flux was maximum or the corresponding change in heat flux for further 
increase or decrease in the plate spacing was not substantial. In fact, the optimum spacing 
is a very important parameter in deciding the thermal design of electronic equipment 
consisting f a number of P(3’s. By thermal design we basically mean arrangement of 
predesignef.' F'rBs in that particular equipment. 

Fig.vre 4.9 shows the variation of optimum spacing with static pressure for 
different p.ate surface temperatures. It can be observed from the figure that for a static 
pressure h.'low 4.0 Pa, the optimum spacing more or less increased with decrease in 
pressure, '.'ht maximum value of the optimum spacing occxuxed for the natural 
convection Ic. when the static pressure was 2 ero. Maximum optimum spacing value was 
45 mm fc- th.e plate surface temperature of 70°C. Below 2.0 Pa, static pressure, the 
optimum snacing increased with increase in plate surface temperature. Above this value 
of static pressure, there was sudden change in trend and the optimum spacing was found 
to decrease v/ith increase in plate surface temperature. This trend of the data seemed to be 
unrealistic and needed further experimentation If the measurement errors are accounted 
for, it may be easy to conclude that the optimum spacing tended to approach some 
limiting value for each plate temperature as the static pressure increased beyond 4.0 Pa. 
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4.6 Effect of Static Pressure on Optimum Heat Flux 


Optimum heat flux refers to heat flux value corresponding to the optimum spacing w^hich 
was the maximum achievable heat flux if the plate spacing was varied along with the 
static prcss’.ue. Figure 4.10 shows the variation in optimum heat flux with change in static 
pressure for different plate surface temperatures. It can be observed that the value of 
optimum heat flux increased with increase in plate surface temperature as well as with 
increase in static pressure. I’or UWT of 50°C, optimum heat flux was found to first 
increase siddimly and then vary linearly with the increase in static pressure. At this 
temperatun; a'^'d static pressure of 15.0 Pa, the maximum value of heat flux was around 
0.17 W/cm\ 

Fo)' th'i: plate surface temperatures of 60°C and l(fC the optimum heat flux 
increased ulmost linearly (starting from 0.19 and 0.29 W/cm^ respectively) with the 
increase in st.'.>ic pressure. Maximum value of heat flux for these two temperatures at 
15.0 Pa was around 0.30 and 0.48 W/cm^ respectively. 

4.7 Correlation 


Suggested correlations for optimum Spacing and optimum heat flux are as follows. 
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15 < X < 70, 50 < T < 70, 0 < P <15 

Xopt is it. m., T is in °C, P is in Pa and Qopt is in W/cm 
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4.8 Error Analysis 

Due to cert,.’ir. inaccuracy in the data acquired, errors in the final results are inevitable. 
Different instruments have different inaccuracies, which get compounded in the final 
results. The Knowledge of inaccuracy of different instruments used in the current 
experiment.il study and the associated errors introduced in the subsequently calculated 
quantities a:e very important before one can really appreciate the experimental results. 
Some of th«.‘ inportant information related to experimental inaccuracies is listed below: 


Table 4.1 


S.No. 

Instruments 

Accuracy 

Associated 

Reading 

Error 

1. 

Microvolt 

j'leter 

+0.5 %fed 

Temperature (“C) 

±0.5 % 

2. 

Vane 

anemometer 

±0.25 %fsd 

Flow velocity (v) 

±0.25 % 

3. 


1.0% 

Plate spacing (m) 

1.0% 

■I 

Temperature 

Controllers 

±2.0 % 

Temperature (“C) 

±2.0 % 


Table 4.2 


S.No. 

Calculated Parameter 

Percent Error 


Nusselt No. 

±1.5% 

2. 

Rayleigh No. 

±5.0 % 

3. 

Heat Flux 

±0.5 % 
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the plate .' pacing indicating some limiting value at higher static pressures. Additional 
cxperirr.en ation can throw more light on this aspect. Optimum Heat flux was found 
to have a higher value for the larger UWT and more or less linear increment with 
increase in .static pressure values. 

4. Two en pineal relations for optimum plate spacing and optimum heat flux were also 
present. -d. 

4,10 ScGp^* for Future Work 

The Unifor.m Wall Temperature condition can be achieved more effectively by dividing 
the heated plate in a larger number of heated strips and thus incorporating a large number 
of heaters. This requires use of a good Data Acquisition and Feedback Control System. 
The unifo.'m plate temperature value requires having a broad range supposedly between 
60°C to loO^'C. The range of plate spacing and the static pressure also need to be 
increased b-yond the range considered. A very good quality work can be pursued if 
different types of temperature profiles are considered over the plate surface. 
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